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Abstract 
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1 Introduction 



The Bethe-Salpeter equation in instantaneous approximation for its interaction kernel 
has proven to represent a very powerful tool for the description of bound states within 
the framework of relativistic quantum field theories. In contrast to the Bethe-Salpeter 
equation, the instantaneous Bethe-Salpeter equation (or "Salpeter equation" ) may be 
formulated as an eigenvalue problem for the "Salpeter amplitudes" which describe the 
bound states under study; its eigenvalues are the masses M of these bound states. In 
the physical sector of bound states of positive norm, all mass eigenvalues M are real 

In a recent series of papers 0, H, ^, we developed a technique for the (approximate) 
determination of the mass eigenvalues M and the corresponding Salpeter amplitudes x 
in an almost entirely analytical way, by conversion of the instantaneous Bethe-Salpeter 
equation into an eigenvalue problem for an explicitly given, analytically known matrix. 
Here we demonstrate, for our (now almost "standard") example, how the ground-state 
mass M for given spin of the bound state may be obtained, with comparable precision, 
even without the necessity of the numerical diagonalization of a (rather large) matrix. 



2 PC = —1 instantaneous Bethe-Salpeter equation 

In order to facilitate an eventual comparison with the results of previous investigations, 
we try to imitate the analysis of Refs. |^ ^ |] to the utmost possible extent. First of all, 
let us base our formalism on the same two simplifying assumptions as in Refs. p|, |^, ^: 

1. Every (full) fermion propagator entering in the Bethe-Salpeter equation may be 
approximated by the corresponding free one, involving a mass parameter which 
then must be interpreted as some effective mass of the bound-state constituents. 

2. The particles forming the bound state under consideration have equal masses m. 

Moreover, let us focus our interest to fermion-antifermion bound states with spin J, 
parity P = (— 1)'^"'"^ and charge-conjugation quantum number C = (—1)'^ (this means 
PC = — 1 for all J), denoted by ^Jj in the usual spectroscopic notation. The Salpeter 
amplitude x describing these states involves two independent components, \E'i and \E'2- 
Its momentum-space representation x(k) is given, for two fermions of equal masses m 
and internal momentum k, in the center-of-momentum frame of the bound state, by 



X(k) 

where 



^i(k)^^^^^ + *2(k)7" 



75 



^(A;) = VA;2 + m2 , A; = |k| , 

denotes the energy of a free particle of mass m and momentum k. Assume that the 
Salpeter amplitude x describes bound states with the total spin J and its projection J3. 
Its two independent components ^i(k) and ^'2(k) may then be factorized, according to 

into the radial wave functions \I'i(A;) and '$2{k), and the spherical harmonics y^mi^) 
for angular momentum i and its projection m; the latter depend on the solid angle Q, 
which encompasses the angular variables, and satisfy the orthonormalization condition 



mm' 
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A crucial point in the construction of any Bethe-Salpeter model for bound states is 
the determination of the Lorentz structure of the Bethe-Salpeter kernel. According to 
the analyses presented in Refs. ^, |^ , for an interaction potential rising linearly with 
the distance of the two bound-state constituents — as is, for instance, frequently used in 
relativistic quark models of hadrons in order to describe the confining quark-antiquark 
interaction arising from quantum chromo dynamics; see, e.g., Refs. [§, |^ — a kernel with 
a pure time-component Lorentz-vector Dirac structure yields stable solutions whereas 
a kernel with a pure Lorentz-scalar Dirac structure certainly does not. For this reason, 
we have chosen in Refs. [|], ^ to discuss pure time-component Lorentz-vector kernels. 

Finally, merely for notational simplicity, we confine ourselves to the case J = 0, i.e., 
to ^So bound states, with the spin-parity-charge conjugation assignment J"^*" = 0""*". 

For pure time-component Lorentz-vector interactions, that is, for a kernel with the 
Dirac structure F^F = 7° 07°, the instantaneous Bethe-Salpeter equation describing 
P C = —1 fermion-antifermion bound states with total spin J = is equivalent to the 
following set of coupled equations for the radial wave functions "^i^k) and '^2{k) ^: 

""dk' k'^ 



-^Voik,k') ^2ik') = M^,ik) 







the interaction between the bound-state constituents, described by the static potential 
V{r) in configuration space, enters in form of the expressions 

00 
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VL{k,k')=8nJdrr^V{r)jL{kr)jL{k'r), L = 0, 1 , (2) 


where jn{z) {n = 0, ±1, ±2, . . .) are the spherical Bessel functions of the first kind [TO 



3 Ground-state mass eigenvalue 

The structure of the PC = —1 instantaneous Bethe-Salpeter equation (|lD suggests to 
solve this set of equations by expressing, for bound-state mass M 7^ 0, from the first of 
Eqs. (0), the Salpeter component \E'i(/c) in terms of the Salpeter component '^2{k), and 
inserting the resulting expression into the second of Eqs. (p. By this procedure the set 
of equations (|l|) is reduced to an eigenvalue equation for \E'2(A;) with as eigenvalue: 



00 



r dk' k'"^ 

M'm^ik) = AE\k)^>2{k) + 2E{k) j^'^^Vo{Kk')^>2{k' 



rdk' k'^ m'V,{k,k') + kk'V,{k.k') 

/w W) 

]dk' k'' m'Vo{k,k') + kk'V,{k,k') J dk"k"^ 

Our goal is to find an as far as possible analytic albeit approximate characterization 
of the bound-state masses M. To this end, we rely on an approximate description of the 
analyzed bound states by (trial) states |0) which involve a real variational parameter 
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Our choice for |0) makes use of the exponential; our |0) is defined in terms of its real 
configuration-space or momentum-space representation 0(r) or respectively, by 

Normahzability of the Hilbert-space states 10) requires /i to be strictly positive: /U > 0. 
These trial functions (/)(r) and are related by the Fourier-Bessel transformations 

<^PP^ 3o{P'r)(t>{p) , <P{p) = y- J drr"^ jo{pr) (f){r) , (4) 



and they satisfy the normalization conditions 

00 00 
J drr^0^(r) = J dpp'^(j)'^{p) = 1 . 


The latter feature allows us to extract, from Eq. (^, the eigenvalue in closed form: 

00 00 00 

= 4 fdkk'^ E\k) (f)'^{k) + - — - / dk E{k) 0(A;) / dA;' k'^ Vo{k, k') (t){k') 





00 , , , r, 00 



00 g 00 

Now, the angular- momentum component Vi(A;, k') of the interaction potential V{r) 
involves — according to its definition (|]) — the spherical Bessel function ji{z) of order 1. 
Consequently, in order to deal with those terms in Eq. (^) which involve the component 
Viik, k'), it is advisable to introduce, in addition, states 1-0) that are related to angular 
momentum 1; the configuration-space and momentum-space representations of are 
then related by a Fourier-Bessel transformation similar to Eq. (^) but involving ji{z): 

Hr) = -^reM-f^r), ^b) = -1 ^^^^ , 

with 

t/j{r)=iJ—J dpp^ ii{j>r)%lj{p) , ip{p) = j drr^ ji{j>r)ip{r) . 



These further functions -(/^(r) and '4){p) are also normalized to unity, that is, they satisfy 

00 

drr^'0^(r) = / dpp^ = 1 • 
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The evaluation of the kinetic energy in the expression (^) for is straightforward: 



dk k' E\k) (j)\k) = dkk^ (j)\k) +m' = fi' + m' . 



However, in order to evaluate the interaction terms in Eq. (|^), we have to approximate 
the various expressions entering in the integrands of these terms in the following way: 



E{k)(/){k) 

^ m 



m 



k(p{k) 

m 



bm , 

cij{k) , 
dilj{k) , 
e(f){k) . 



Expressed in terms of the coefficients b, c, d, and e defined above, the square (|^) of the 
bound-state mass eigenvalue M of the instantaneous Bethe-Salpeter equation (|l]) may 
then be traced back to the expectation values 



= {<i)\V{r)\(P) = J drr^V{r)(l)\r) , 



oo 

= {'ilj\V{r)\^lj) = [ drr^V{r)'ilj^{r) 



of the interaction potential V{r) with respect to the approximation functions (j) and i/j 



(6) 



For the factors b, c, d, and e analytic expressions may be found; the latter will obviously 
involve the mass m of the bound-state constituents and the variational parameter /i: 



dkk^E{k) 4>\k) =^x^ 
dk k^ 



1 X 



(2;2- 1)5/2 l^^^S 

3 - 16 ^2 + 48 



d 



2^3 

^/3 

dkk^^*{k)(P{k) = i — /i , 



[x^ 



1)9/2 



2 x^ 



ISvrxio 



[k) = 



X — Ax^ + 







(3.44; 









with the Gauss hypergeometric series F, defined, in terms of the gamma function F, by 

m 



FM 



T{u)Tiv) 



E 

n=0 



Tiu + n) Tiv + n) z' 



Yiw + n) 
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and the abbreviation 

m 

X = — 

for the ratio of the mass parameters involved. It is rather easy to convince oneself that 
these expressions for the coefficients 6, c, d, and e reduce, in the (ultrarelativistic) limit 
m — > 0, to the lowest entries of the corresponding matrices given in Ref. 0, where the 
instantaneous Bethe-Salpeter equation has been studied for the simpler special case of 
vanishing masses of all bound-state constituents, and, for the particular value n = m, 
to the lowest entries of the corresponding (similarly defined) matrices given in Ref. 0, 
where, because of the nonvanishing masses of the bound-state constituents, the desired 
analytical treatment forced us to fix the value of the variational parameter /i to /i = m. 
The necessarily numerical optimization of the right-hand side of the analytic result 
with respect to the variational parameter /i then yields a first approximation to the 
masses M of bound states described by the instantaneous Bethe-Salpeter equation (|l|). 



4 Accuracy of the approximation 

The decisive question clearly is whether it is possible to achieve a satisfactory accuracy 
of the bound-state mass M calculated from the approximation represented by Eq. (^). 
The answer to this question will, of course, depend on the interaction potential V{r). 
The general expressions for the expectation values of V{r) for arbitrary power-law 



potentials may be deduced from Refs. 12, O. Here we study as a simple but, from 



the physical point of view, nevertheless relevant example the case of a linear potential: 

V{r) = \r, A>0. 

Inserting the required expectation values V^^^ and V^^^ of this potential (which may be 
found, e.g., in Refs. @, H), 

T/(0) = — T/(l) = — 

into our general formula (||) for the mass squared of the lowest ^Sq bound state yields 

= A{m' + fi') + — {b + m-" e) +—C* d + ^ (s + 5 \c\') . (7) 

Minimization of this expression with respect to n, for m = 0.1 GeV and A = 0.2 GeV^, 
gives M = 1.703 GeV; this is only 2.5% larger than the "exact" value M = 1.661 GeV 
computed in Ref. ^ by the diagonalization of a 25 x 25 matrix (cf. Table 1 of Ref. [0). 
Table |] compares the squared masses predicted by Eq. (|^) with the findings presented 
in Table 1 of Ref. Again, the relative errors are only of the order of a few percent. 



5 Summary and conclusion 

The demand to derive, also for the case of massive bound-state constituents, analytical 
matrix representations of the instantaneous Bethe-Salpeter equation led us, in Ref. 0, 
to identify the variational parameter in the trial states used in Ref. with the mass m 
of the bound-state constituents. However, the rate of convergence of the corresponding 
bound-state masses M with increasing matrix size is, because of this identification, for 
small but nonvanishing masses m of the bound-state constituents not extremely rapid. 
For a first idea of the location of M, one might want to avoid matrix diagonalizations. 
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Table 1: Lowest mass eigenvalues M of the instantaneous Betlie-Salpeter equation, 
with a time-component Lorentz- vector interaction kernel, describing bound states with 
spin-parity-charge conjugation assignment J^^ = ^ of two spin-| fermions of mass 
m, which experience a confining interaction described by a linear potential V{r) = Xr 
with slope A = 0.29 GeV^ The values of arising, for various values of the mass m 
of the bound-state constituents, within the present approach (third column) differ by a 
few percent (fourth column) from the values of obtained in Ref. by expanding 
the Salpeter amplitude x over a basis of 25 harmonic-oscillator eigenfunctions (second 
column). Moreover, the relative errors decrease with increasing constituents' mass m. 



m [MeV] 


M2 [GeV2] 


M2 [GeV2] 


Relative Error 


Relative Error 




(Ref. i) 


(present work) 


of M2 [%] 


of M [%] 


300 


4.357 


4.568 


4.8 


2.4 


500 


5.255 


5.478 


4.2 


2.1 


900 


8.248 


8.515 


3.2 


1.6 



Here, by applying the same variational technique as used in Ref. , we have been able 
to demonstrate that, at least for ground states, the bound-state mass M may be found 
with reasonable accuracy, certainly sufficient for the initial steps of a fitting procedure. 
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